Detections of C2H, cyclic-CgHa, and H^^CN in NGC 1068 
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ABSTRACT 

We used the Nobeyama 45-m telescope to conduct a spectral line survey in the 3-nim band 
(85.1-98.4 GHz) toward one of the nearest galaxies with active galactic nucleus NGC 1068 and 
the prototypical starburst galaxy NGC 253. The beam size of this telescope is ^ 18", which was 
sufficient to spatially separate the nuclear molecular emission from the emission of the circumnu- 
clear starburst region in NGC 1068. We detected rotational transitions of C2H, cyclic-C3H2, and 
H-'^'^CN in NGC 1068. These are detections of carbon-chain and carbon-ring molecules in NGC 
1068. In addition, the C2H N = 1-0 lines were detected in NGC 253. The column densities of 
C2H were determined to be 3.4 x lO^^ cm'^ in NGC 1068 and 1.8 x 10^^ cm'^ in NGC 253. The 
column densities of cyclic-C3H2 were determined to be 1.7 x 10^'^ cm~^ in NGC 1068 and 4.4 x 
10^"^ cm^^ in NGC 253. We calculated the abundances of these molecules relative to CS for both 
NGC 1068 and NGC 253, and found that there were no significant differences in the abundances 
between the two galaxies. This result suggests that the basic carbon-containing molecules are 
either insusceptible to AGN, or are tracing cold (Trot ^ 10 K) molecular gas rather than X-ray 
irradiated hot gas. 
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Introduction 



To date, about 40 molecular species have been 
identified in nearby external galaxies (e.g. list 
available in CDM^il; Miiller et al. 2005). As a re- 
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suit, it is possible to study molecular abundances 
and chemical reactions in external galaxies. For 
example, a difference in molecular abundance has 
been reported between the nearby starburst galax- 
ies NGC 253 and M 82 (e.g. Mauersberger & 
Henkel 1991, Takano et al. 1995), and the rea- 
son for this difference has been discussed. They 
suggested that the temperature in NGC 253 is 
higher than that in M 82 and/or that the size of 
the high-temperature region in NGC 253 is much 
larger than that in M 82. In addition, Takano 
et al. (2002) pointed out the peculiarity of the 
molecular abundance in M 82 among nearby star- 
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burst galaxies. Some groups have indicated that 
the specificity of the molecular abundances of M 82 
is because of the chemistry of photon-dominated 
regions (PDRs) (e.g. Garcia-Burillo et al. 2002; 
Martin et al. 2009). 

Different galaxies have different properties and 
activities. These include those with small or large 
amounts of gas, having a starburst and/or ac- 
tive galactic nucleus (AGN), and with interac- 
tion, etc. Molecular line observations of these 
different galaxies allow us to study the effects of 
these different properties/activities on the molec- 
ular medium. In fact, some groups have sug- 
gested that it is possible to diagnose power sources 
in dusty galaxies using molecular line ratios (e.g. 
Kohno et al. 2001, Uscro et al. 2004, Kohno 2005, 
Imanishi et al. 2007, Kohno et al. 2008, Krips et 
al. 2008). The observation of the molecular gas 
chemistry of the AGN toward NGC 1068, one of 
the nearest galaxies with an AGN, has been re- 
ported by Usero et al. (2004), Perez-Beaupuits 
et al. (2009), and Garcfa-Burillo et al. (2010). 
They observed molecular lines of HCO, H^^CO+, 
SiO, HCO+, HOC+, and CN, and concluded that 
the circumnuclear disk of NGC 1068 is a giant X- 
ray dominated region (XDR). There are theoret- 
ical studies on the XDR model by Lcpp & Dal- 
garno (1996), Maloney et al. (1996), and Mei- 
jerink et al. (2005; 2007). However, further sys- 
tematic observations of molecular lines are indis- 
pensable to study the impact of AGN on the inter- 
stellar medium. Nevertheless, no systematic unbi- 
ased scans of millimeter /submillimeter molecular 
lines toward AGN have been published to date, 
although such scans do exist for two nearby star- 
burst galaxies (NGC 253, Martin et al. 2006; M 
82, Naylor et al. 2010). Therefore, we started 
a project to conduct a line survey in the 3-mm 
band of NGC 1068 using the 45-m telescope at 
Nobeyama Radio Observatory (NRO) in Japan. 
The beam size of this telescope (18" at 86 GHz) 
is smaller than the size of the circumnuclear star- 
burst ring in NGC 1068 (d - 30"; e.g. Planesas et 
al. 1991) , and it is therefore essential to study the 
impact of the AGN on the surrounding molecules; 
this will enable us to mitigate the contamination 
of the molecular lines from the circumnuclear star- 
burst region in NGC 1068. We also observed NGC 
253 to compare the effects of AGN on molecular 
abundance. This is an on-going project, and we 



present here the initial results of the survey. 

In this study, we report detections of the ro- 
tational transition lines of C2H, cyclic-C3H2, and 
Hi-'^CN in NGC 1068, and detection of the C2H 
N = 1-0 transition in NGC 253 as well as other 
previously reported lines. 

The C2H molecule was first detected in inter- 
stellar clouds by Tucker et al. (1974). However, 
extragalactic detections are limited to just a few 
objects (M82, Henkel et al. 1988; NGC 4945, 
Henkel et al. 1990; NGC 253, Martin et al. 2006). 
The C2H molecule is important to study the for- 
mation and the characteristics of carbon-chain 
molecules, because it is related to the carbon- 
chain growth. The cyclic-C3H2 molecule has been 
observed in a large number of galactic sources. 
In external galaxies, it has been detected in Cen 
A (Seaquist & Bell 1986; Beh & Seaquist 1988), 
NGC 253, and M82 (Mauersberger et al. 1991). 
Information on the abundances of basic carbon- 
containing molecules is important to understand 
the carbon chemistry in NGC 1068. The detec- 
tion of H-'^^CN in the external galaxies has only 
been reported toward NGC 253 by Mauersberger 
& Henkel (1991). It is interesting that the isotope 
of HCN was detected in NGC 1068, because the 
isotope ratio gives the constraint on the optical 
depth of the HCN line; this will help us to unveil 
the nature of the over-luminous HCN emission at 
the center of NGC 1068. 

2. Observations 

The observations at the 3-mm region (85.1 - 
98.4 GHz) toward NGC 1068 and NGC 253 were 
carried out from February to May, 2009 and from 
January to May, 2010 using the 45-m telescope. 
The total observational times under good weather 
condition are 44.5 and 22.5 hours for NGC 1068 
and NGC 253, respectively. The adopted central 
position and the systemic velocity (Vlsr) of NGC 
1068 are aj2ooo = 02 h 42 m 40.798 s, ^72000 = -00° 
00' 47.938" and 1150 km s-\ respectively (Schin- 
nerer et al. 2000). Similarly, the adopted central 
position and Vlsr of NGC 253 are 0172000 = 00 h 
47 m 33.3 s, Sj2aoo = -25° 17' 23 " (Martin et al. 
2006) and 230 km s'^ (Takano et al. 1995), re- 
spectively. In particular, we focused on the cyclic- 
C3H2 Jk^.k^ = 2i.2 lo.i {i^rcst = 85.338906 GHz), 
Hi^CN J= 1-0 (i/rest = 86.340167 GHz, F=2-l), 
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CzH TV = 1 (t'rcst = 87.316925 GHz, J = 3/2- 
1/2, F = 2-1), HCN J = 1-0 (t-rest = 88.631847 
GHz, F = 2-1), CH3OH Jk = 2k-Ik {vr^^t = 
96.741377 GHz, Jk = 2o Iq A+), and CS J= 2-1 
(I'rest = 97.980953 GHz) transitions. These rest 
frequencies were taken from the database of Lovas 
(2004). 

The dual-polarization sideband-separating re- 
ceiver (Nakajima et al. 2008) was used with typi- 
cal system noise temperatures including an atmo- 
sphere of about 150-300 K, and the image rejec- 
tion ratio (IRR) was typically better than 10 dB. 
It was checked by the IRR measurement system 
(Nakajima et al. 2010), which employs an artifi- 
cial signal injection from the top of the receiver op- 
tics. The half-power beam widths (HPBWs) at 86 
GHz were 18.3" and 18.4" for H-polarization (H- 
pol.) and V-polarization (V-pol.), respectively, in 
2009, and 19.2" and 18.9" for H-pol. and V-poL, 
respectively, in 2010. The measured main beam ef- 
ficiencies (r/mb) were 0.43 ± 0.02 and 0.42 ± 0.02 
for H-pol. and V-pol., respectively, in 2009, and 
0.49 ± 0.04 and 0.47 ± 0.03 for H-pol. and V-pol., 
respectively, in 2010. The backend used consisted 
of eight digital spectrometers (Sorai et al. 2000), 
each having a bandwidth of 512 MHz and a res- 
olution of 605 kHz. We used four spectrometers 
for each polarization, lined up in the direction of 
frequency. 

The line intensities were calibrated by the 
chopper wheel method. The telescope pointing 
was checked by observing the nearby SiO maser 
sources, o-Cet for NGC 1068, and R-Aqr for NGC 
253 every ^-^1.5 hour. The pointing accuracy was 
better than about 5". The position-switching 
mode was employed for the observations where 
the reference position was taken at the position 
with an azimuthal angle difference of +5' for both 
galaxies. The integration time on each position 
was 20 sec. 

3. Results 

The main results include detections of C2H, 
cyclic-CsHa, and H^^CN toward NGC 1068 and 
the C2H N = 1 transition toward NGC 253. In 
addition, we observed the HCN J = 1-0 line and 
marginally detected CH3OH {Jk = 2^-1 k) to- 
ward NGC 1068. Very recently, Garcia-Burillo et 
al. (2010) reported the detection of the CH3OH 



group of transitions {Jk — 2^ Ik and Jk — 3j(— 
2k)- Moreover, detection of CS {J = 2-1) has 
already been reported by Tacconi et al. (1997). 
Nonetheless, we detected this transition in NGC 
1068 using a single dish telescope. Figures 1 and 2 
show the observed line profiles, and Table 1 sum- 
marizes the parameters derived from the Gaussian 
fits to the observed lines. 

3.1. C2H 

The C2H = 1-0 line consists of 6 hyperfine 
components (J = 3/2-1/2, = 1-1, F = 2-1, F 
= 1-0 and J = 1/2-1/2, F = 1-1, F = 0-1, F = 
1-0). We detected two fine structure components 
because of the line blending. Thus, the fitting has 
been carried out with a comb of Gaussian profiles 
at the rest frequencies of the hyperfine components 
(Gottlieb et al. 1983) using the same widths and 
with fixed optical depths (r) ratios relative to the 
main component based on the theoretical relative 
intensities (Tucker et al. 1974). The width of the 
components and the optical depth were taken as 
free parameters. 

The main-beam brightness temperature (Tmb) 
can be calculated from the following equation: 



Tn^b = r]u{J{Te.) - J{Tbg)}{l 



where 



1 



^ ^ k exp(|f 



(1) 



(2) 



1) 

The beam filling factor, = / (^s + ^b)' ac- 
counts for the dilution eff'ect due to the coupling 
between the source and the telescope beam in the 
approximation of a Gaussian source distribution of 
size 9s (FWHM) that is observed with a Gaussian 
beam size 0b (HPBW). The adopted 6^ of NGC 
1068 is 4", which is assumed from the distribution 
of HCN (Heifer & Britz 1995; Krips et al. 2008). 
The adopted 9^ of NGC 253 is 20", which is as- 
sumed from the distribution of CS (Martin et al. 
2006). We observed the C2H lines only in 2009, 
and 0h was measured to be ~ 18.4" at 86 GHz in 
2009. Thus, for the C2H observation, the values 
of r]b{ for NGC 1068 and NGC 253 with the 45- 
m telescope are 0.045 and 0.54, respectively. We 
assumed the excitation temperatures (Tex) toward 
NGC 1068 to be 7.6 K, which is derived from the 
analysis of the rotation diagram of the CS lines 
(see subsection 3.3), and that toward NGC 253 
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to be 5.8 K, which is derived from the analysis 
of the rotation diagram of the C2H hnes (see af- 
ter the next paragraph). The employed cosmic 
background radiation (Tbg) is 2.7 K. The fitting 
results were overlaid on the spectra (Figure 2 (a) 
and (b)). The Gaussian profiles for C2H in NGC 
253 do not reproduce the observed spectrum well. 
The peak position of the main component seems to 
be shifted possibly because of the self-absorption 
effect. Because of this shift, the overall fitting re- 
sults are not as good as those of NGC 1068. 

We obtained the optical depth of the C2H main 
component {N = 1 0, J = 3/2 1/2, F = 2 1) of 
NGC 1068 and NGC 253 to be 0.06 and 0.05, re- 
spectively. The total of the optical depth of all 
the C2H components in both galaxies are 0.15 and 
0.13, respectively. This indicates that the C2H N 
= 1-0 lines have small (<C 1) optical depth in NGC 
1068 and NGC 253. If the source size of NGC 1068 
is changed to 10", to be able to see the effect of the 
assumed source size, then Tex becomes 8.5 K. In 
this case, the optical depth is calculated to be 0.01. 
Therefore, the result of the small optical depth of 
C2H is rather robust considering the uncertainty 
of the assumed source size. The obtained column 
densities of C2H in NGC 1068 and NGC 253 are 
3.4 X 10^^ cm~^ and 1.8 x 10^^ cm~^, respec- 
tively. For this calciilation, we adopted a dipole 
moment of 0.8 debye (Tucker et al. 1974). On 
the other hand, the column density of 1.2 x 10^^ 
cm~2 was reported in NGC 253 from the C2H N = 
2-1 line by Martin et al. (2006). Thus, our result 
compares well to theirs. 

Because the C2H N = 2-1 line was observed 
by Martin et al. (2006), we made a rotational di- 
agram of C2H toward NGC 253 (Figure 3). We 
adopted the basic equation for the rotation di- 
agram (see for example Tiuner 1991). The fit 
to the iV = 1 (this work) and 2 1 (Martin et 
al. 2006) transitions of C2H toward NGC 253 ap- 
pears to indicate a low excitation gas of 5.8 
K. This value is close to the T^t of NO (6 K) 
and NS (7 K), which are reported by Martin et 
al. (2006). These molecular lines trace gas of a 
low-temperature component. 

3.2. cyclic-C3H2 

Spectra of cyclic-C3H2 are shown in Figure 2 (c) 
and (d) for NGC 1068 and NGC 253, respectively. 
The 2i^2-lo,i transition was clearly detected in 



both galaxies. Figure 3 shows the rotational di- 
agram toward NGC 253, where we employed the 
dipole moment of 3.43 debye reported by Kanata 
et al. (1987) (Brown et al. (1987) reported it to 
be 3.32 debye). The statistical weight is 3 for the 
ortho levels and 1 for the para levels due to two 
equivalent hydrogen nuclei. In this diagram, we 
assumed an ortho-to-para ratio of 3. The fit to 
the 2i^2-lo,i (ortho; this work) and 32,2-2i_i (para; 
Martin et al. 2006) transitions of cyclic-C3H2 to- 
ward NGC 253 appears to indicate a low excitation 
gas with Trot of 7.6t J j K. Martin et al. (2006) 
detected three transitions of the cyclic-C3H2 lines 
and obtained Trot of 9 ± 8 K from their rotational 
diagram. However, they reported that the lines 
are blended with other lines except the 32,2-2i_i 
transition. This is the reason why we used only 
their 32,2-2i^i transition in our rotation diagram. 
As a result, a more probable Trot was obtained. 
The column density of cyclic-C3H2 toward NGC 
253 is derived using the value of the intersection 
of the linear regression with the y-axis of Figure 3 
and the partition function of cyclic-C3H2 at 7.6 K. 
For the calculation of the partition function, the 
energy levels in Vrtileck et al. (1987) were used. 
The obtained column density is 4.4 x 10^"^ cm~^. 
For NGC 1068, the obtained column densities are 
1.7 x 101=^ cm-2 (6)3 = 4") and 1.8 x 10" cm'^ 
{0s = 10"). These column densities in NGC 1068 
were derived from the above analysis of the ro- 
tation diagram assuming Tot of 7.6 K, which is 
derived from the rotation diagram of CS (see the 
next subsection). 

3.3. CS 

Figure 1 (d) shows the CS J= 2-1 spectrum to- 
ward NGC 1068. In addition, the CS J= 3-2, 5-4, 
and 7-6 transitions have already been reported by 
Bayet et al. (2009) in NGC 1068. Using these 4 
lines and assuming 9s of 4", the rotational temper- 
ature was calculated to be 12.6 K from the rotation 
diagram as inducted by a dashed line (Figure 3). 
However, Bayet et al. mentioned that the J= 7-6 
line is marginally detected. Subsequently, with the 
exception of the J = 7-6 line, the performed fits 
are shown by the thin solid line in Figure 3, and 
the obtained Tot is 7.6 K. The result of the fitting, 
excluding the J = 7-6 line, is deemed appropriate, 
because the fitting error is significantly small. If 
the 0s of NGC 1068 is changed to 10", then the 
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Trot will be 13.0 K and 8.5 K (excluding the J 
= 7-6 line). We find that the slope of the linear 
regression in the CS rotational diagram with our 
new data point is still consistent with previous re- 
sults (Bayet et al. 2009), indicating that physical 
properties of the gas traced by CS (J = 2-1) are 
similar to those traced by the higher J transitions 
of CS. Trot of 7.1 K (excluding the J = 7-6 line) 
has been reported by Bayet et al. (2009), and our 
result is consistent with their result. 

3.4. HCN and H^^CN 

The relative intensities of the H^^CN (Figurel 

(a)) and HCN (Figurel (b)) lines can be used as a 
measure of the optical depth of HCN. The optical 
depth is derived using the following equation, as- 
suming the same excitation temperature for HCN 
and Hi^CN, and a i^c/^C isotopic ratio of 50 
(Lucas & Liszt 1998): 

TX(Hi3CN) _ 1 - e-^(HCN)/50 
TX(HCN) ~ i_e-^(HCN) • (3) 

We obtained an optical depth of about 2.6 for the 

HCN J = 1-0 line toward NGC 1068. Thus, the 
HCN J = 1-0 hue in NGC 1068 is optically thick. 

4. Discussion 

It is interesting to examine the effects of AGN 
on molecular abundance. For example, is there a 
significant difference in the molecular abundance 
in circumnuclcar disk because of the presence of an 
AGN? Here, we focused on the C2H and cyclic- 
C3H2molecules that were detected in this study 
toward NGC 1068. The relative abundances of 
C2H were calculated with respect to CS in both 
galaxies to study the effects of the AGN. CS was 
employed to calculate the relative abundances be- 
cause the multi-transition data of CS lines have 
been published (Bayet et al. 2009). In addition, 
the critical density of CS is about 10"*"^ cm~^, and 
it traces the region where the densities are higher 
than that traced by CO. The column densities of 
CS toward NGC 1068 and NGC 253 were obtained 
as reported in subsection 3.3. 

For C2H, we obtained the abundances relative 
to CS of 5.5 (6*3 = 4") and 5.7 (6*^ = 10") in NGC 
1068, and 4.3 in NGC 253. For cyclic-C3H2, we ob- 
tained the abundances relative to CS of 0.03 {6s = 
4") and 0.16 {6^ = 10") in NGC 1068, and of 0.10 



in NGC 253. Martin etal. (2006) reported the col- 
umn density in NGC 253 to be 3.0 x lO^^ cm'^, 
and a relative abundance calculated to be 0.07 in 
the same way as shown above. It is almost equal 
in ratio for both of the molecules, and we found no 
significant differences in the relative abundances of 
NGC 1068 and NGC 253. Thus, there is a possi- 
bility that these basic carbon-containing molecules 
are insusceptible to AGN. Another possibility is 
that C2H and cyclic-C3H2 mainly exist in an ex- 
tended gas away from AGN, because the rotation 
temperature of C2H (7.6 K) is lower than the ki- 
netic temperature of the gas around the central 
few hundred parsecs (> 70 K; Tacconi et al. 1994). 
However, as far as we know, there is no theoretical 
study on C2H and cyclic-C3H2 in XDR. Therefore, 
such model calculations are essential. 

5. Conclusions 

We started line surveys in the 3-mm band to- 
ward NGC 1068 and NGC 253. The main re- 
sults include detections of C2H, cyclic-C3H2, and 
H13CN toward NGC 1068 and the C2H N = 1- 
transition toward NGC 253. We calculated the 
relative abundances of C2H and cyclic-C3H2 with 
respect to CS in both galaxies and found no signif- 
icant differences in the results between NGC 1068 
and NGC 253. Thus, it is concluded that these 
basic carbon-containing molecules are insuscepti- 
ble to AGN and/or these molecules exist in a cold 
gas away from AGN. 

We continue an unbiased line survey toward 
NGC 1068 and NGC 253 and examine the abun- 
dances of other molecules in order to study the ef- 
fects of AGN on the interstellar medium. A clear 
separation of the circumnuclcar disk from the star- 
burst ring is important to understand the abun- 
dances of the molecular gas around AGN. With 
the advent of ALMA, we will be able to study the 
AGN chemistry with much higher spatial resolu- 
tion. 

This work is based on observations with the 
NRO 45-m telescope, as a part of the line sur- 
veys of the legacy projects. The authors thank 
the project members for helpful discussion. We 
also thank all of the staff of the 45-m telescope 
for their support. We arc also grateful to Ryohei 
Kawabe for his advice and support to the line sur- 
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vey project. 

Note : After the submission of this manuscript, 
we found two relevant papers. Snell et aL de- 
tected C2H in NGC 1068 (astro-ph 1101.0132). It 
was vc;ry recently piiblished in The Astronomical 
Journal 141 38. Costagliola et al. detected C2H 
and H13CN in NGC 1068 (astro-ph 1101.2122). 
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Fig. 1.— (a) H^^CN J= 1-0, (b) HCN J= 1-0, 
(c) CH3OH Jk = 2k-Ik, and (d) CS J= 2-1 lines 
toward NGC 1068 obtained with a velocity resolu- 
tion of 20 km s~^. The first order baselines were 
subtracted. The reference frequencies are listed 
in the section of observations. The gray curves 
overlaid on the profiles are the results of Gaussian 
least-squares fits. Note that the CH3OH line is at 
best tentatively detected (see text). 




NGC 10 


38 c-C,H, 


... 1 ... 1 . 


(c) 



400 800 1200 1600 

V,,Jkms-'] 



NGC 253 



A. J 



3 



K,Jkms-'] 



(b) 

'466 800 

c'-C,H,' 



(d)| 
too 



Fig. 2.— The C2H N = 1-0 lines toward (a) 
NGC 1068 and (b) NGC 253, and the cychc-CgHa 
JKa,,Ka = 2i^2^1o,i lines toward (c) NGC 1068 and 
(d) NGC 253. For details, see the caption of Fig- 
ure 1. 
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E / k[K] 

Fig. 3. — Rotational diagrams of C2H, cyclic- 
C3H2 in NGC 253, and CS in NGC 1068. The 
black and open circle symbols represent the results 
of C2H of this work and the N = 2-1 transition 
by Martin et al. (2006), respectively. The black 
and open triangle symbols represent the results of 
cyclic-C3H2 of this work and the 32,2~2i^i transi- 
tion by Martin et al. (2006), respectively. The 
black and open square symbols represent the re- 
sults of CS of this work and the J = 3-2, 5-4, and 
7-6 transitions by Bayet et al. (2009). The fit 
with all the data of CS is shown by the dashed 
line, and the fit excluding the data point of the J 
= 7-6 line, is shown by the thin solid line. 
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Table 1: The line parameters 



line 


transition 




^mb 


^SR 


At;(FWHM) 


rms (Tmb) 






[mK km s~^] 


[mK] 


[km s-i] 


[km s-i] 


[mK] 


NGC 1068 














cyclic-C3H2 


JKa,Kc = 2i,2-lo,i 


1200 


5 


1148 


248 


1.7 




J= 1-0 


1900 


9 


1111 


234 


4.1 


C2H 


N= 1-0 


6900 


19 


1139 


209 


3.3 


HCN 


J= 1-0 


35000 


128 


1117 


336 


19.1 


CH3OH 


Jk = '^K-^K 


1500 


7 


1104 


235 


3.8 


CS 


J= 2-1 


6400 


30 


1174 


312 


4.8 


NGC 253 














cyclic-C3H2 


jRa.Ka — 2l,2-l0,l 


5500 


30 


272 


243 


4.6 


C2H 


N= 1-0 


48000 


129 


219 


199 


9.1 
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